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Abstract 

After extending the Clarkson-Kruskal's direct similarity reduction ansatz to a more gen- 
O^l ' eral form, one may obtain various new types of reduction equations. Especially, some 

r~ — ^ lower dimensional turbulent systems or chaotic systems may be obtained from the general 

i^ \ type of similarity reductions of a higher dimensional Lax integrable model. Especially, the 

t^^ , Kuramoto-Sivashinsky equation and an arbitrary three order quasi-linear equation which 

O ■ 

includes the Korteweg de-Vries Burgers equation and the general Lorenz equation as two 

^^ , special cases are obtained from the reductions of the (2-t-l)-dimensional dispersive long wave 

^ ■ equation system. Some types of periodic and chaotic solutions of the (2-|-l)-dimensional 

Ph ■ dispersive long wave equation system are also discussed. 

> 

X. 

c3 ■ 1 Introduction 

To reduce a higher dimensional nonlinear physical model to some lower dimensional ones is 
one of the most important approaches in the study of nonlinear science. Usually one use the 
standard Lie group approach to reduce a higher dimensional partial differential equation (PDE) 
to lower dimensional ones[||]. Lately, the so-called nonclassical Lie group analysis is established 
to find lower dimensional similarity reductions [^. To find some lower dimensional reductions by 
using the classical and nonclassical Lie group approaches, one has to use some tedious algebraic 
procedures. In the past decade, to avoid the tedious algebraic calculation in the finding of 
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the similarity reductions, a simple direct powerful method is developed|3|, Q. Using the direct 
method, various new similarity reductions of many physical models are found though these 
reductions can also be obtained lately from the nonclassical Lie group approach p, y, 0]. In Q, 
the direct method is extended to find some types of conditional similarity reductions which have 
not yet been obtained by means of the present classical Lie group approach and nonclassical 
Lie group approach. In this paper, we try to extend the direct method in another direction to 
find lower dimensional reductions that may not be obtained by using the present classical and 
nonclassical Lie group approaches. 

In the next section, we discuss the general aspect on the direct reduction method. In section 

3, the (2+l)-dimensional dispersive long wave equation (DLWE) is used as a concrete example 
to realize new reduction idea and to find some new lower dimensional reductions. In section 

4, we use some numerical solutions of the lower dimensional reduction models to discuss some 
types of exact solutions of the (2+l)-dimensional DLWE. The last section is a short summary 
and discussion. 

2 General reduction ansatz of direct method 

To reduce many types of (n+l)-dimensional nonlinear PDEs, 

A(xj, u, Ux^, UxiXj, •••, i,j = 0,1, ■■■in) = A[u\=0, {xo = t), (1) 

it is proven that the special ansatz 

U = a{xo,Xi,...,Xn) +P{xo,Xi,...,Xn)w{Co,^l, ■■., Cn-l), (2) 

is sufficient instead of 

U = U{xo,Xi, ...,Xn,w{Co,Cl, •••, ^n-l)), (3) 

where w{S,j) = w{S,o,(,i, ..., ^n-i) satisfies an n-dimensional PDE. In (2) and (3), u and w may 
be some multi-component fields. However, to reduce a higher dimensional PDE to some lower 
dimensional ones one may use some more general ansatzs instead of (3). For instance, the ansatz 
(3) may be extended as 

u = Uixi, w^^, w^^^^^^, ..., u'c,^^,^...^,^ = U[w]. (4) 

In other words, some types of derivatives of the reduction function may be included in the 
primary reduction ansatz. However, to find some concrete results is quite difficult by using the 



general ansatz (4). By using the similar procedure of the simplification from (3) to (2) for many 
types of significant mathematical physics models, one may simplify (4) to 

u = Uq [w] + Ui [w] w^^^ (.^^ . . .^^.^ , (5) 

where w^. ^. ,,,^. is one of the highest derivatives of w included in (4) while C/o[t(;] and [/i[w] 
are wc. c^ p. independent. 

3 Special new reductions of the (2+l)-diniensional DLWE 

To give out some concrete results from above general reduction ansatz, we take the (2+1)- 
dimensional dispersive long wave equation (2DDLWE) 

Uyt + VXX + UxUy + UU^y = 0, (6) 

rit + Ux+ V^x + u-q^ + Uxxy = (7) 

as a simple example. The equation system (6) and (7) is first obtained by Boiti et al. Q as 
a compatibility condition for a 'weak' Lax pair. The infinite dimensional Kac-Moody-Virasoro 
type symmetry structure of the model is revealed by Paquin and Winternitz|p!o|. The more 



general Woo symmetry is given in|ll|. It is proven that[|l2[ the 2DDLWE system is fails in 
passing the Painleve test both at the WTC's (Weiss- Tabor-Carnevale) |jl^ meaning and at the 
ARS's (Ablowitz-Ramani-Segur) meaning[0]. Using the special ansatz (2), nine types of two 
dimensional similarity reductions and thirteen types of ODE (ordinary differential equation) 
reductions has been given by one of the present authors (Lou)[p^]. 

For simplicity further, we taking the reduction ansatz (5) in a specific form, 

U = Fi{t,y,w)Wx + Fo{w,Wx,Wxx) + F2{w,Wx)Wxxx, (8) 

V = r] + l = Uy = Fiy{t,y,w)wx, ■w = w{x,t). (9) 

The reason why we take the ansatz (8) is that we try to find the reduction equations have the 
following three order autonomous PDE form 

Wt = aWxxx + F3{w, Wx, Wxx) (10) 

for some possible functions F3. The ansatz (9) degenerates two equations (6) and (7) to a same 
one. 

Substituting (8)- (10) into (6) and/or (7) yields 

Fiy{t,y,w){F2{w,Wx)wx + a)wxxxx + f{t,y,w,Wx,Wxx,Wxxx) = 0, (11) 



where f{t,y,w,Wx,Wxx,Wxxx) = / is a complicated expression of the indicated variables. Be- 
cause of / is Wxxxx independent, (11) is valid only for 

F2{w, Wx) = -awx^. (12) 

Substituting (12) into (11), we have 

Fiy{t,y,w){F3i^^^{w,Wx,Wxx) + 1 + WxFoitj^^{w,Wx,Wxx))Wxxx + fl{t,y,W,Wx,Wxx) = 0, (13) 

where fi{t,y,w,Wx,Wxx) = /i is independent of Wxxx- From eq. (13) we immediately have 

Fsiw, Wx, Wxx) = -Wxx-WxFo{w, Wx, Wxx) + F^iw, Wx). (14) 

By using Eq. (14), (13) is simplified further to 

{2wxFiy{t,y,w)Fi{t,y,w) + Fiy{t,y,p)F3u,Aw,Wx) + 2wxFiyyj{t,y,p))wxx 

+f2{t,y,w,Wx) = 0, (15) 

where f2i't,y,w,Wx) = /2 is Wxx independent. Integrating (15) once with respect to y, we have 

(wxFi{t,y,wf + Fi{t,y,w)F3^^{w,Wx) + 2wxFi^{t,y,w) + h{t,w,Wx)j Wxx 

+f2{t,y,w,Wx) = 0, (16) 

with f3{t,w,Wx) being an integrating function. Because of the Wx independence of Fi{t,y,w), 
by vanishing the first term of (16), we get 

F3{w,Wx) = F32{w)wl + F3o{w), fi{t,W,Wx)=WxFi{t,w), (17) 

and 

2Fi^{t,y,w) + Fi{t,y,wf + 2F32{w)Fi{t,y,w) + Fi{t,w) = {). (18) 

Because of (17) and (18), Eq. (16) is simplified finally to 

Fu{t,y,w) + Fi{t,y,w)F3Q^{w) - ^Fso{w){F,{t,y,w) + 2F32{w))Fi{t,y,w) + F5{t,w) = 0.(19) 
The compatibility condition of (18) and (19) requires that 

F5{t,w) = -F30^u,{w) + F32{w)F30u,H - -F4{t,w)Fso{w) + F32U'w)F3o{w) (20) 



and 



F4{t,w)F30yj{w) - F32{w)F3o{w)F32wiw) - 2F30w{w)F32wiw) + -F3o{w)F4i,{t,w) 
-F32wwiw)F3o{w) - F32{wfF30i,{w) + -Fsott,^^^, (ttj) + -F4t{t,w) = 0. (21) 



Now the final results show us that the 2DDLWE (6) and (7) possesses the following reduction 

Wt = aWxxx - Wxx - WxFo{w, Wx, Wxx) + F32{w)wl + F3o{w) (22) 

with arbitrary functions Fq{w, Wx, Wxx), -p32(u^) and F3q{'w) and 

u = -aw~'^Wxxx + Fo{w, Wx, Wxx)+ Fi{t,y,w)wx, (23) 

V = r]+l = Fiy{t,y,w)wx, (24) 

where Fi{t, y, w) is determined by two compatible Riccati equations (18) and (19) while -F4(t, w) 
and F^{t, w) are determined by (20) and (21). The simplest solution of (18)-(21) reads 

F4{t,w) = Fs2{w) = F5{t,w) + 2A2 = 0, (25) 

Fsoiw) = A2w'^ + Aiw + Ao, (26) 

Fi{t,y,w)= ^ (27) 

w + q{y,t) 

qt{y, t) = Aiq{y, t) - Ao - A2q{y, tf, (28) 

where Aq, Ai, A2 are arbitrary constants. 

Prom the reduction equation (24), we can see that though the original 2DDLWE system 
is Lax integrable, possesses infinitely many symmetries and abundant multi-soliton structures, 
there still exist various nonintegrable lower dimensional reductions because of the entrance 
of three arbitrary functions Fq{w, Wx, Wxx), -f32(^i') and F30. For instance, if we select 
Fo{w, Wx, Wxx), F32iw) and F3o{w) simply as 

- WxFo{w, Wx, Wxx) + F32{w)wl + F3o{w) = WWx, (29) 

then (22) becomes the well known KdV-Burgers equation 

Wt = aWxxx - Wxx + WWx (30) 

which is one of the possible candidate to describe the turbulence phenomena in fluid physics and 
plasma physics [p^ , |T7| . If the functions Fq{w, Wx, Wxx), F32{w) and F30 are fixed to satisfy 

-WxFow,^{w, Wx, Wxx) + F32{w)wl + F3o{w) 

= — [wxxWx + (c + l)wl] - w'^Wx - (6 + c)wxx " wc{b -ba + w"^) (31) 

w 

with Q = —1 and a, b, c are arbitrary constants, then (24) becomes a (l+l)-dimensional 
extension 

Wt = -Wxxx - {b + c + l)wxx H [wxxWx + (c + l)wl] - w'^Wx - wc{b -ba + w"^) (32) 

w 



of the famous chaotic system, the Lorenz system |18| 

Ws = -c{w - g), gs = {a- h)w - g, hg = wg - bh. (33) 

Actually, the travelling wave reduction of (32), w = w{x + b{c+l)t) = w{s), is totally equivalent 
to the Lorenz system (33). 

In principle, any order of derivatives of w may be included in the ansatz (5). And some types 
of more complicated reduction equations can be obtained. For instance, if we insert a fourth 
order derivative Wxxxx term into the reduction ansatz, we may obtain many fourth order (1+1)- 
dimensional PDE reductions. Here we list only a special example for the reduction equation has 
a famous Kuramoto-Sivashinsky (KS) equation form||n 



2wx{ai + asQ) (03 =F ^)wxx , 

u = ± — h aiw 

oo + aiw + a2Q + a^^wQ Wx 

H {a^Wxxxx + (aiC2 + a2,CQ)u? + (c2ao + cioi + a2Co + 02)^ + ciOq), (34) 

Wx 

_ QyPxiasao -Q2ai) f„^. 

ao + aip + a2Q + aspQy 

Qt = aoco + (ciai + a2Co — C2ao)Q + (cios — C2a2)Q^, (36) 

wt + aiwwx + a2W + a^Wxx + aiWxxxx = 0, (37) 

where Oj, Oi, ? = 0, 1, ...,5 are arbitrary constants and cq, ci and C3 are arbitrary functions of 
t. Various interesting properties of the chaotic KS equation (37) have been studied by many 
authors, say, |19| and the references therein. 



In the reduction results (22) and (37), the independent variables are simply taken as x 
and t. Actually, extending these independent variables to some more general forms is possible 



because the model possesses infinitely many symmetries with some arbitrary functions [^, 11 1. 



For instance, using the finite transformations given in |1C], all the independent variables of the 



systems (22) and (37) are changed to some general forms naturally. 

4 Special solutions 

Now an interesting question is which kinds of exact solutions can be obtained from our new 
reduction equations? In this section, we write and plot down some interesting exact solutions. 

4.1 Multi-dromion solutions 

If we take Fq{'w, Wx, Wxx) as 

Fo{w, Wx, Wxx) =Qw - w~^{wxx - wIF32{w) - F3o{w)), (38) 




Figure 1: Four dromion solution of the 2DDLWE for the field w = 77 + 1 given by (24) with (27) and (41)-(44). 
then we know that the w equation (22) is just the well known KdV equation 



wt = awxxx — (JwWd 



(39) 



Then we can use the N soliton solutions of the (l+l)-diniensional KdV equation to construct 
the multi-dromion solutions by taking 



N 



q{y, t) =ao + Y^ a„ tanh(/„y - Vn) 



(40) 



n=l 



with Oj, J = 0, 1, ..., N being arbitrary constants and Aq = Ai = A2 = 0. Fig. 1 is a plot of the 
four dromion solution with w being two soliton solution of the KdV equation 



w = -2a(ln(l))xx, 

(p = 1 + ex.p{kix + akit + xi) + ex.p{k2X + 0^2^ + X2) 

+71 r^J exp((A;i + k2)x + a{kf + k^)t + xi + X2) 

(fci + k2y 



and 



(41) 



(42) 



q = 4 + tanh(^iy - yi) + 2tanh(Z2y - 2/2) 



(43) 



while the other constants are fixed as 



ki = 1, ^2 = 1-1, ^1 = ^2 = 1, a = -1, xi = -3, X2 = 3, yi = -3, 2/2 = 3. 
The figure (1) and all other figures of this paper are plotted at time t = 0. 



(44) 



4.2 Periodic and chaotic line soliton solutions 

If q is still given by (43) while w is given by (33), then we may obtain some kinds of periodic or 
chaotic line soliton solutions. Because no one has given ever out any exact explicit solutions of 
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Figure 2: (a) Periodic two solution of the Lorenz system, (b) The related periodic two line soliton solution of 
the 2DDLWE for the field v. 

(33) , we can only use the numerical solutions of the generalized Lorenz system to construct exact 
solutions of the 2DDLWE. For several types of parameter ranges, the solutions of the Lorenz 
system are periodic while for other types of parameter ranges, the solutions of the Lorenz system 
are chaotic ones. Fig. 2 is a plot of the periodic two line soliton solution of the 2DDLWE with 
the parameters of the Lorenz system (33) is fixed as 



and 



a = 350, 6= -, c= 10 
3' 



gr = 200 + tanh y 



(45) 



(46) 



From Fig. 2b we can see that the line solution is localized in y direction and periodic in 
s{= X + b{c + l)t) direction when the parameters are selected appropriately as (45). 

Fig. 3 plots the chaotic line soliton solution of the 2DDLWE with the parameters of the 
Lorenz system (33) given by 



a = 60, b= -, c= 10 



(47) 



while the q function still given by (46). Obviously, Fig. 3 shows us that when the parameters 
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Figure 3: (a) chaotic solution of the Lorenz system, (b) The related chaotic line soliton solution of the 2DDLWE 
for the field v. 

of (33) are located at the chaotic regions, the corresponding solution becomes a chaotic straight 
line soliton solution which is localized in y direction and chaotic in s direction. 



4.3 Space periodic and chaotic solutions 

From (28) we know that in some cases, the function q may be an arbitrary function of y, so 
we may also select it as a solutions of the Lorenz system (33) with the replacement of the 
independent variable s ^ y. When the function q = q{y) and w = w{s) = w{x + 6(c + l)t) are 
all the solutions of the Lorenz system, we can obtain many types of solutions which are periodic 
or chaotic in both directions. 

Fig. 4 is a plot of a periodic solution of the 2DDLWE which has periodic property in both 
directions while q{y) and w{s) are chosen as both the solutions of the Lorenz system (33) with 
(45). 

Fig. 5 is a plot of an exact solution of the 2DDLWE which is periodic in y direction and 
chaotic in x direction. The corresponding solutions for q and w are all determined by (33) but 
with different parameters (45) and (47) respectively. 

Fig. 6 shows a chaotic solution of the 2DDLWE in both directions. The related solutions for 
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Figure 4: Periodic two solution of the 2DDLWE in both direction for the field v. 




Figure 5: Plot of the exact solution of the 2DDLWE for the field which is periodic in y direction and chaotic in 
X direction. 
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Figure 6: A typical solution of the 2DDLWE which is chaotic in both directions. 
q and w are all determined by (33) with the same parameters (47). 

5 Summary and discussions 

In summary, the CK's direct similarity reduction ansatz is extended to a much more general form. 
Using the general reduction ansatz, one may obtain various new lower dimensional reduction 
equations including many turbulence and chaotic systems. Taking the 2DDLWE as a concrete 
example, and a slightly special reduction ansatz with three order derivatives of the reduction 
field, we obtain a general three order quasi-linear equation, which includes the KdV, MKdV, 
KdV-Burgers and the generalized Lorenz system as special examples, as a special reduction 
of the 2DDLWE. The known KS system and other types of higher order models may also be 
obtained from the reductions of the 2DDLWE. 

The reductions (30), (32) and (37) are known as some typical turbulence and chaotic systems 
while the 2DDLWE is known as an 1ST integrable model. The reason why some lower dimen- 
sional turbulence and chaotic systems can be reduced from a higher dimensional integrable 
(under some particular meanings) model is that for a higher dimensional integrable model, some 
types of lower dimensional arbitrary functions do enter into its general solution. 

Using the solutions of the lower dimensional models we may obtain many kinds of new 
solutions for the 2DDLWE. Especially, using the numerical solutions of the Lorenz systems, 
some types of periodic line solitons, chaotic line solitons, periodic-periodic solution, periodic- 
chaotic and chaotic-chaotic solutions of the 2DDLWE can be obtained. 



In Ref. [20 1, using the variable separation approach, we have also pointed out that the 
turbulence and chaotic systems can be obtained from other "integrable" models like the Davey- 
Stewartson equations and the asymmetric Nizhnik-Novikov-Veselov equation because of the 
entrance of the lower dimensional arbitrary functions in the general solutions. The more about 

11 



the method and the effects of the reduced turbulence system on the original model(s) are worthy 
of studying further. 
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Huang and C-p Sun. The work was supported by the National Outstanding Youth Foundation 
of China, the Research Fund for the Doctoral Program of Higher Education of China and the 
Natural Science Foundation of Zhejiang Province of China. 
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